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Abstract. We propose a method of optimally controlling the tradeoff of speed 
and fidelity of state transfer through a noisy quantum channel (spin-chain). This 
process is treated as qubit state-transfer through a fermionic bath. We show that 
dynamical modulation of the boundary-qubits levels can ensure state transfer with the 
best tradeoff of speed and fidelity. This is achievable by dynamically optimizing the 
transmission spectrum of the channel. The resulting optimal control is robust against 
both static and fluctuating noise in the channel’s spin-spin couplings. It may also 
facilitate transfer in the presence of diagonal disorder (on site energy noise) in the 
channel. 
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One dimensional (ID) chains of spin-| systems with nearest-neighbor couplings, 
nicknamed spin chains, constitute a paradigmatic quantum many-body system of the 
Ising type [T]. As such, spin chains are well suited for studying the transition from 
quantum to classical transport and from mobility to localization of excitations as a 
function of disorder and temperature [2], In the context of quantum information (QI), 
spin chains are envisioned to form reliable quantum channels for QI transmission between 
nodes (or blocks) piSJ- Contenders for the realization of high-fidelity QI transmission 
are spin chains comprised of superconducting qubits 0E], cold atoms [3 El El EDJ, 
nuclear spins in liquid- or solid-state NMR [TTJ [12, d3| [H, T5] 1T6], |T7, [18], quantum dots 
m, ion traps |20. [21] and nitrogen-vacancy (NV) centers in diamond [221 '23, [231 M] • 

The distribution of coupling strengths between the spins that form the quantum 
channel, determines the state transfer-fidelities [31 [Mi 2?. j2Hi Ml I3U] . Perfect state- 
transfer (PST) channels can be obtained by precisely engineering each of those couplings 
[3T1 [32l 271 [28l 129 , 1551 1 301 [53] , Such engineering is however highly challenging at 
present, being an unfeasible task for long channels that possess a large number of control 
parameters and are increasingly sensitive to imperfections as the number of spins grows 
[35, |31l ESI S3- A much simpler control may involve only the boundary (source and 
target) qubits that are connected via the channel. Recently, it has been shown that 
if the boundary qubits are weakly-coupled to a uniform (homogeneous) channel (i.e., 
one with identical couplings), quantum states can be transmitted with arbitrarily high 
fidelity at the expense of increasing the transfer time j38j ]39l 40; 0H 12. 33) @31 36]. Yet 
such slowdown of the transfer may be detrimental because of omnipresent decoherence. 

To overcome this problem, we here propose a hitherto unexplored approach for 
optimizing the tradeoff between fidelity and speed of state-transfer in quantum channels. 
This approach employs temporal modulation of the couplings between the boundary 
qubits and the rest of the channel. This kind of control has been considered before for a 
different purpose, namely to implement an effective optimal encoding of the state to be 
transferred [33] • Instead, we treat this modulation as dynamical control of the boundary 
system which is coupled to a fermionic bath that is treated as a source of noise. The goal 
of our modulation is to realize an optimal spectral filter [361 sa SHI 1321 EH EH E21 EH E3] 
that blocks transfer via those channel eigenmodes that are responsible for noise-induced 
leakage of the QI [55] . We show that under optimal modulation, the fidelity and the 
speed of transfer can be improved by several orders of magnitude , and the fastest possible 
transfer is achievable (for a given fidelity). 

Our approach allows to reduce the complexity of a large system to that of a 
simple and small open system where it is possible to apply well developed tools of 
quantum control to optimize state transfer with few universal control requirements 
on the source and target qubits. In this picture, the complexity of the channel is 
simply embodied by correlation functions in such a way that we obtain a universal, 
simple, analytical expression for the optimal modulation. While in this article we 
optimize the tradeoff between speed and fidelity so as to avoid decoherence as much 
as possible, this description [361, 37, IS E- EE El, Si EH Ml 133.. allows one to actively 
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suppress decoherence and dissipation in a simple manner, since it may be viewed as a 
generalization of dynamical decoupling protocols [ 561 [57 , 581 I59 J. In what follows, we 
explicitly deal with a spin-chain quantum channel, but point out that our control may 
be applicable to a broad variety of other quantum channels. 


1. Quantum channel and state transfer fidelity 


1.1. Hamiltonian and boundary control 

We consider a chain of IV+ 2 spin-| particles with XX interactions between nearest 
neighbors, which is a candidate for a variety of state-transfer protocols MElEnillM 
unum H3UHUI5UI6UI2UI8UIZUI9U20U2I1E21E3U2H [251 [261 EZl EH [291 [30U311EU [321 [33] . 
The Hamiltonian is given by 


H = H 0 + H bc (t), 


Hr) = 


N -1 

X 


J% { 

0 V 


°i a i+1 


+ erf a. 


f+i) , H bc (t) = aft) Y 


Jr 


a i a i+1 


+ ° V i CfV i+ 1 ) , 


( 1 ) 

( 2 ) 


i =1 i£{0,N} 

where H 0 and H bc stand for the chain and boundary-coupling Hamiltonians, respectively, 
af' yi are the appropriate Pauli matrices and J* are the corresponding exchange- 
interaction couplings. 


1.2. Mapping to a few-body open-quantum system 

The magnetization-conserving Hamiltonian H can be mapped onto a non-interacting 
fermionic Hamiltonian |60| that has the particle-conserving form 

n-i 

H o=Y^ (4c l+ i + dc\ + 1 ) , H hc {t) = ait) Y J (4d+i + c*4+i) , (3) 

i =1 ie{0,N} 

where Cj = \e 1 ^ rJ ^ (Ji crj create a fermion at site j and a ± = a x ± ia v . The 
Hamiltonian H 0 can be diagonalized as H 0 = Yk=i u kb\i>ki where b\ = (j\uk)(f 
populates a single-particle fermionic eigenstate \u>k) of energy oj k , and | j) = |0..01 J 0..0) 
denote the single-excitation subspace. Under the assumption of mirror symmetry of the 
couplings with respect to the source and target qubits Jj = Jn-i, the energies 04 are 
not degenerate, c <4 < 04 + 1 , and the eigenvectors have a definite parity that alternates 
as 04 increases [28]. This property implies that ( 7 ( 04 ) = (—l) fc-1 (IV — j + 1 ( 04 ) and 
allows us to rewrite the boundary-coupling Hamiltonian as 

N N 

H bc (t ) = a(t) J 0 cl'Y,(M UJ k)bk + ot(t)J N c^ N+1 'Yf—l) k ~ 1 (N\ujk)bk + h.c. (4) 

k= 1 k= 1 

For an odd N, there exists a single non-degenerate, zero-energy fermionic mode in 
the quantum channel, labelled by k = z = jjjgj [ 25 ^ A s a consequence, the 
two boundary qubits (0 and IV + 1) are resonantly coupled to this mode. Therefore, 
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we consider these three resonant fermionic inodes as the “system” S and reinterpret the 
other fermionic modes as a “bath” B. In this picture, the system-bath SB interaction 
is off-resonant. Then, we rewrite the total Hamiltonian as 

H = H s {t) + H B + H SB {t), (5) 


where 

N 

HB = ^2 Ukblbk, Hs{t) = s + (t)j z b z + h.c., (6) 

k^z,k=l 

HsB(t)=s+(t) ^2 Jkbk + S-(t) *22 Jkbk + h.c., ( 7 ) 

k£.kodd even 

with s±(t) = a(t)(c' 0 ± cjv +1 ), J k = Ji(l|o; fc ), k odd = {1,3, provided k odd ^ z, and 

keven = {2, 4, IV — 1}. 

The form (|5]) is amenable to the application of optimal dynamical control of the 
multipartite system (3H1 El, E21 EB, E3J : such control would be a generalization of 
the single-qubit dynamical control by modulation of the qubit levels [48, HH 150, 51] E2, 
1531, -Tfj . To this end, we rewrite Eq. (JT]) in the interaction picture as a sum of tensor 
products between system Sj and bath Bj operators (see Appendix A) 


= < 8 ) 

3 = 1 

From this form one can derive the system density matrix of the system, ps(t), in the 
interaction picture, under the assumption of weak system-bath interaction, to second 
order in Hsb, as jl6] [48J 

6 

Ps(t ) = Ps( 0) - t ^2 h'Psi 0)] + h-c., (9) 

i,i '=1 


Appendix A 


where 


1 4 pt 

RiAt) = 7 J2 dt> 
1 . - /0 




( 10 ) 


with dv y ( 7 ) = Tr s \^Bj{r)B-i (O)pe(O)} denoting the correlation functions of bath 
operators and Hj,j(f) being a rotation-matrix in a chosen basis of operators hi used to 


represent the evolving system operators, Sj(t) = 'YfVtj^(t)t)i (Appendix A). The solution 


2—1 


(J9|) will be used to calculate and optimize the state-transfer 


fidelity in what follows. 


1.3. Fidelity derivation 

We are interested in transferring a qubit state |^o) initially stored on the 0 qubit to 
the N + 1 qubit . Here j^o) is an arbitrary normalized superposition of the spin-down 
|0 0 ) and spin-up |1 0 ) (single-spin) states. To assess the state transfer over time T, we 
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Figure 1. (Color online) (a) Top: State transfer through a spin-channel with 
boundary-controlled couplings. Bottom: Boundary-controlled spin chain mapped to a 
non-interacting spinless fermionic system (dashed rectangle) which couples to the bath 
fermionic-modes k (red even k and blue odd k lines) with strengths The two 

boundary spins 0 and N +1 are resonantly coupled to the chain by the fermionic-mode 
^ with a coupling strength J z a(t) (green lines), (b) Bottom: Interacting bath-spectrum 
G±(oj) including noise effects (grey color) bounded by the Wigner-semicircle (maximal- 
disorder) lineshape (dashed contour) with a central gap around uj z . In the central gap, 
an optimal spectral-filter I 7 t,±(w) is shown (green color). Top: generated 

by boundary-control a p (t ): p = 0 (black dotted), p = 1 (violet dashed), p = 2 (orange 
thin). The red vertical lines are the nearest bath-spectrum eigenenergies. The inset 
is a zoom on the tails of the filter spectrum that protects the state transfer against a 
general noisy bath with a central gap. 

calculate the averaged fidelity F(T) = _|_ F.n+i(T) _|_ i jgj w hich is the state- 

transfer fidelity averaged over all possible input states l^o)- In the interaction picture, 
fo,N+i(T) = \ s (i>\ps{T ) |-0) s | where | ip) s = |lo) IO^Oat+i) and |^) s <g) \ip) B is the 
initial state of S + B. 

In the ideal regime of an isolated 3-level system, perfect state transfer occurs when 
the accumulated phase due to the modulation control 

<t>{T) = J z [ a(t)dt (11) 

Jo 

satisfies <f>(T) = Obviously, this condition does not strictly hold when the system- 
bath interaction is accounted for, yet it is still adequate within the second-order 
approximation in Hsb used in Eq. (|9]) . In this approximation, fo,N+i(T) takes the 
form 


fo,N+i(T) = 1-CCH, 


(12) 
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Figure 2. (Color online) Infidelity 1 — F(T) as a function of transfer time T under 
optimal control a p (t): p = 0 (black dotted), p = 1 (violet dashed), p = 2 (orange thin), 
when the noise or bath-correlation spectrum G[u>) is bounded by the Wigner-semicircle 
with a central gap around io z (Appendix C ). 


where 


C (t) = 5ft f dt [ dt?y'n±(t)n±(if)$±(t -t)). 

Jo Jo ± 


(13) 


Here, $±(£) = V\ c , I JJ 2 e l0Jkt are the bath-correlation functions, while £l+(t) = 

v 7 L — J KtzK 0 dd(even) v 7 

a(t)cos(y/2(j)(i)) and f l-(t) = aft) are the corresponding dynamical-control functions 
(Appendix A Appendix B). In the calculations we considered \iP)b — |0)_b- However, 


in the weak-coupling regime the transfer fidelity remains the same for a completely 
unpolarized state (651 HU or any other initial state of the bath. 


In the energy domain, Eq. (13) has the convolutionless form | 481 H9l [50, [51] [52, 5^ 


/ OO 

dwJ2F T ,±(u) G ±(u), 

-OO _L 


where the Fourier-transforms 


G±(u) = — dt$±(t)e luJ \ F Tj ±(w) = — | / dtto±(t)e 


icot 12 


2 vr 


2tt ' 


(14) 


(15) 


are the bath-correlation spectra, G±(oj), associated with odd(even) parity modes and 
the spectral filter functions, Ft,±(oj), which can be designed by the modulation control. 


2. Optimization method 

To ensure the best possible state-transfer fidelity, we use modulation as a tool to 


minimize the infidelity £(T) in (13 14) by rendering the overlap between the interacting 
bath- and filter-spectrum functions as small as possible mm- 
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The minimization of C(T) in (13) can be done for a specific bath-correlation function 


of a given channel which represents a non-Markovian bath. The Euler-Lagrange (E-L) 
equation for minimizing ( (T) with the energy constraint 


E(T) = Jf \a(t)\ 2 dt 
Jo 


turns out to be 


- A - <k - A = o 

dt [ d<j> dp dp ’ 


(16) 


(17) 


where A is the Lagrange multiplier and 0 = J z a. The optimal modulation can be 
obtained by solving the integro-differential equation 


4>{t) = — 


pEQ(t,<f>(t),f(t)) 


Jz\JIf d t\Jo | 2 ’ 


(18) 


where 


= fo dt '®(t - t')&U ( d * + a t,] cos(V2(/)(t))cos(V2(l)(t')) + ^ d f,) 


+ |^($ + (O)co S 2 (v / 20(t)) + <h_(O)). 

(19) 

The solution of Eq. (18) should satisfy the boundary conditions 0(0) = 0 and 0(T) = 
to ensure the requireastate transfer. 

In general the bath-correlations have recurrences and time fluctuations due to 
mesoscopic revivals in finite-length channels. Therefore, it is not trivial to solve Eqs. 


(18 19) analytically rather than solving them numerically for each specific channel. We 


however are interested in obtaining universal analytical solutions for state-transfer in the 
presence of non-Markovian noise sources. To this end, we here discuss suitable criteria 
for optimizing the state transfer in such cases. 

We require the channel to be symmetric with respect to the source and target 
qubits and the number of eigenvalues to be odd. These requirements allow for a central 
eigenvalue that is invariant under noise on the couplings. This holds provided a gap 
exists between the central eigenvalue and the adjacent ones, i.e. they are not strongly 
blurred (mixed) by the noise, so as not to make them overlap. At the same time, we 
assume that the discreteness of the bath spectrum of the quantum channel is smoothed 
out by the noise, since it tends to affect more strongly the higher frequencies [34l [36 . 37j. 
Then, if we consider the central eigenvalue as part of the system, a common characteristic 
of G±(cu) is to have a central gap (as exemplified in Fig. 0 1 )- 

Therefore, in order to minimize the overlap between G-t(u;) and Ft.±(uj) for general 


gapped baths, and thereby the transfer infidelity in (14), we will design a narrow 
bandpass filter centered on the gap. 

We present a universal approach that allows us to obtain analytical solutions for 
a narrow bandpass filter around u z . Since G_( oj) has a narrower gap than G + (ui), 
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we optimize the filter Ft~(uj) under the variational E-L method. We seek a narrow 
bandpass filter, whose form on time-domain via Fourier-transform decays as slowly as 
possible, so as to filter out the higher frequencies. This amounts to maximizing 


Ft At) = 


F T -(to) e lUT du} = / a{t)a{t + r)dt, 


( 20 ) 


subject to the variational E-L equation (17), upon replacing ( by F T _. Since there is no 
explicit dependence on 0, the second term therein is null, Ft - — A eE) = 0, yielding 


a(t + r) + a(t - r) = A£«(£) + A^, 


( 21 ) 


where A^ is the Lagrange multiplier and A^ is an integration constant chosen to satisfy 
the boundary conditions obeyed by the accumulated phase 0- 

Analytical solutions of (21) are obtainable for small r, corresponding to the 
differential equation 

b(t) = -A E Oi{t) + A (22) 


with A e — 


-Qe~ 2 ) 


and A«£ = -4 • It has a general solution 


a (t) = Asin(u v t ) + Bcos(uj v t ) + C. 


(23) 


The unknown parameters are then optimized under chosen constraints, e.g. on the 
boundary coupling, the transfer time, the energy, etc. 

The frequencies uj v that give a low and flat filter Ft A A outside a small range 
around u> = oj z = 0 are oo v = ?reZ, since the components of a(t) that oscillate with 
oj v then interfere destructively. Only if n = 0,1,2 will the filter have a single central 
peak around cu = 0, and the contribution of larger frequencies will be suppressed, while 
the filter-overlap with the central energy level will be maximized; for larger n, the central 
peak splits and additional peaks appear at larger frequencies. 

Therefore, the analytical expressions for the optimal solutions satisfying 0(0) = 0 
and 0(T) = ^ are found to be 


where p — 0,1, 2, 


a p (t) = aMsin p 



7r 


OlM Cp 


V2 J Z T 


(24) 

(25) 


and Cp — TiL. l ap — 1, Ci — |, C 2 — 2). Here p — 0 means static control, while 
p — 1,2 stand for dynamical control. Note that T and olm = max{a p (t)} cannot 
be independently chosen. If the transfer time is fixed, then the maximum amplitude 
depends on p, ce M = according to Eq. (25). Similarly, if the maximum amplitude 
is kept constant, then the transfer time will depend on p, T = T p . by Eq. (25). 
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The different solutions in Eq. (24) are sinc-like bandpass filter functions around 

0 that become narrower as T increases. For p = 0, which satisfies the minimal-energy 

2 

condition E min (T Q ) = the corresponding filter is the narrowest around 0, but it has 
many wiggles on the filter tails (Fig. §>) which overlap with bath-energies that hamper 
the transfer. In contrast, the p — 1,2 bandpass filters are wider (for the same T) and 
require more energy, E\ = ^ E min and E 2 = § E min respectively, but these filters are 
flatter and lower throughout the bath-energy domain. 

Hence, the bandpass filter width {i.e. full width at half maximum) and the overlap 
of its tail-wiggles with bath-energies as a function of T, determine which modulations 
a p (t ) are optimal, as shown in the inset of Fig. 0> (F Tt+ (uj) filters out a similar spectral 
range). The shorter T, the lower is p that yields the highest fidelity, because the central 
peak of the filter that produces the dominant overlap with the bath spectrum is then 
the narrowest. However, as T increases, larger p will give rise to higher fidelity, because 
now the tails of the filter make the dominant contribution to the overlap. As shown in 
Fig. [2j the filter for p = 1,2 can improve the transfer fidelity by orders of magnitude in 
a noisy gapped bath bounded by the Wigner-semicircle, which is representative of fully 
randomized channels 


(Appendix C). 


2.2. Optimizing the modulation control for a Markovian Bath 

We next consider the worst-case scenario of a Markovian bath, where the bath- 
correlation functions < h±('r) vanish for r > 0. This is the case when the gap is closed 
by a noise causing the bath energy levels to fluctuate faster than the system dynamics. 
We note that, finding optimal solutions for the noise spectrum of a Markovian bath is 
important for the case where the gap is reduced or even lost in static cases. 


The infidelity function (13) that must be minimized when the correlation time 
Tc = 0, i.e. <&±(t) = S(r), is 

" T _ 

(26) 


ccn = 9f if dt^M($ + ( 0 )cos 2 (V 2 <l>(t)) + MO)). 

Jo Jz 


is now 


The E-L equation under energy constraint © , i 

4>(t) ^<h + (O)cos 2 (v^20(t)) + < l ) _(O)— 2XJ 2 j —\/2 4> 2 {t) ^ + (0)cos(V2(f>(t))sin(V2(j>(t)) =0. 

(27) 

This equation has a non-trivial analytical solution and the modulation that minimizes 
C(T) is given by the following transcendental equation 

T f 0 Ht) yj cos(2v / 20) < f>+(O) + M0) + 2 M0) - 2A J 2 z d(j) 

-t J^ T) ^/2(Mo mMmTmmM)# = 0. 

The infidelity for this optimal modulation almost coincides with the one obtained for 


(28) 


static control (a p= o(t) =«« from Eq. (24)), i.e. 


1 - /••('/') 


n 2 N 


;( 1 - 


7 T 2 N 


6V2JT 16\/2 JT 


T = 


7T 


y/N 


2oi mJ 


(29) 
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Figure 3. (Color online) Transfer infidelity 1 — F(T) for a modulated boundary- 
controlled coupling a p (t) = Qj}SOT p (f) as a function of (a) the transfer time T, (b) 
the maximum value of the boundary coupling p = 0 (empty circles), p = 1 (filled 
squares), p = 2 control (empty diamonds). The quantum channel is a homogeneous 
spin-chain with TV + 2 = 31 spins and J = 1. 


and they only differ by about 0.1%. This optimal modulation can be phenomenologically 
approximated by 

.. . f 7T,, f)l , . 

ait) ~ aajvf + bsin q (—)), q ~ 3.5, - ~ a ~ 0.84, (30) 

T a 3 

assuming no constraints (A = 0). An example of the performance of this solution is 
discussed below and shown in Fig. [5j 


3. Optimal control of transfer in a homogeneous spin-chain channel 


Consider a uniform (homogeneous) spin-chain channel, i.e. Ji = J in Eq. ([Tj) , whose 
energy eigenvalues are u>k = 2Jcos(-Ayy) [58] . In Fig. 3j we show the performance of 


the general optimal solutions (24) for this specific channel as a function of and T. 


The approach based on Eq. 

[a M < i) im m ssi ns eqi m, eu. 


(|13|) strictly holds in the weak-coupling regime 
In this regime (marked with arrows in Fig. 


we found that the transfer time is T, n 


T ydV 


and the infidelity decreases by 


^ ‘P ~ P 2a M J 1 

reducing aM according to a power law, aside from the oscillations due to the discrete 

nature of the bath-spectrum (see Appendix C). The filter tails are sinc-like functions, so 


that when a zero of the filter matches a bath-energy eigenvalue, the infidelity exhibits 
a dip. Aside from oscillations, the best tradeoff between speed and fidelity within this 
regime is given by the optimal modulation with p = 2 (for the system described in Fig. 

§>)■ 


However, this approach can also be extended to strong couplings aM, since it 
becomes compatible with the weak-coupling regime under the optimal filtering process 
that increases the state fidelity in the interaction picture [5HE21I63]. The bandpass filter 
width increases as T decreases; consequently, in the strong coupling regime ~ 1) 


















Optimized dynamical control of state transfer through noisy spin chains 


11 


the filter may now overlap the bath energies closest to uj z , but still block the higher 
bath energies, which are the most detrimental for the state transfer [Ml [36 , 157] . Then, 
the participation of the closest bath energies yields a transfer time T p ^c p jj. There is 
a clear minimal infidelity value at the point that we denote as a° A p / p which depends on p 
(Fig. [§)); thus extending the previous static-control (p = 0) results, where an optimal 
ol°m was found [26] E2, DEJ [36]. The infidelity dip corresponds to a better filtering-out 
(suppression) of the higher energies, retaining only those that correspond to an almost 
equidistant spectrum of Uk around oj z , which allow for coherent transfer [36]. 

Figure [3}o shows that by fixing aw, the dynamical control (p = 1,2) of the 
boundary-couplings reduces the transfer infidelity by orders of magnitude only at the 
expense of slowing down the transfer time T p at most by a factor of 2, ^ ^ < 2. If 

the constraint on can be relaxed, i.e. more energy can be used, the advantages of 
dynamical control can be even more appreciated for both infidelity decrease and transfer¬ 
time reduction by orders of magnitude, as shown in Fig. [3^i. Hence, our main result 
is that the speed-fidelity tradeoff can be drastically improved under optimal dynamical 
control. 


4. Robustness against different noises 


We now explicitly consider the effects of optimal control on noise affecting the coupling 
strengths, also called off-diagonal noise, causing: Jj —> Ji + JjAj(t), i = 1 , ...,1V with 
Aj being a uniformly distributed random variable in the interval [—£j,£j\. Here £j > 0 
characterizes the noise or disorder strength. When Aj is time-independent, it is called 
static noise , as was considered in other state-transfer protocols [Ml EDI SUES]. When 
A i(t) is time-dependent, we call it fluctuating noise [7T]. These kinds of noises will 
affect the bath energy levels, while the central energy oj z remains invariant [MET]. In 
the following we analyse the performance of the control solutions obtained in Sec. [2] 


for these types of noise and later on, in Sec. 4A we discuss briefly the effects of other 
sources of noise. 


J h l. Static noise 

Static control on the boundary-couplings can suppress static noise [34~ . 36] but here we 
show that dynamical boundary-control makes the channel even more robust, because 
it filters out the bath-energies that damage the transfer. To illustrate this point, we 
compare the effect of modulations a p {t) with ccm = a°f^ for p = 0 and 2 in the strong¬ 
coupling regime (Fig. &)• There is an evident advantage of dynamical control with 
p = 2 compared to static control (p — 0), at the expense of increasing the transfer time 
by only a factor of 2, [A ~ 2. In the weak-coupling regime, if we choose such that 
the transfer fidelity is similar for p = 0 and p = 2, then both cases are similarly robust 
under static disorder, but the modulated case p = 2 is an order of magnitude faster. 
Remarkably, because of disorder-induced localization [72] ESI EH d, regardless of how 
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Figure 4. (Color online) Transfer infidelity for a modulated boundary-controlled 
coupling a p (t) as a function of the perturbation strength Sj of the noisy homogeneous 
channel, averaged over N av noise realizations for (a) static and fluctuating noisy 
channels = 0.6, a 0 ff 2 = 0.7 and N av = 10 3 . In static noisy channels, the 
infidelity obtained under static control p = 0 (empty circles) is shown to be strongly 
reduced when dynamical p = 2 control is applied (empty squares). A fluctuating noisy 
channel is less damaging; in the Markovian limit, where the correlation time of the 
noise fluctuations t c —>• 0 (p = 0, green solid circles), the infidelity converges to its 
unperturbed value. The homogeneous channel has N + 2 = 31 spins and J = 1. 
(b) Same plot for fluctuating noise, ranging between static and Markovian noise for 
«m 0 = 0.1, T = 88 and N av = 200. Here r c is the correlation time of the noise 
fluctuations (see text). Faster fluctuations reduce the noise effect and thereby the 
fidelity decay. The inset shows the scaling of the infidelity, where the effective noise 
strength is scaled with y = yj 2r c (all curves overlap). The scaled noise strength 
depends on the noise correlation time r c . 


small is ctM, the averaged fidelity under static noise cannot be improved beyond the 
bound 


1 — F oc Ne 2 j, (ej -C 1). 


(31) 


4-2. Markovian noise 


The worst scenario for quantum state transfer is the absence of an energy gap around 
oj z . This case corresponds to Markovian noise characterized by (Aj(t)Aj(t + r)) = 
S(t), where the brackets denote the noise ensemble average, or equivalent to a bath 
correlation-function that vanishes at r > 0. In this case there is an analytical solution 
for the optimal modulation given by Eq. (28), although the infidelity achieved by it 


almost coincides with the one obtained by the static (p = 0) optimal control (Fig. [5]). 
Counterintuitively, arbitrarily high fidelities can be achieved for such noise by decreasing 
max\a(t)\ and thereby slowing down the transfer. This comes about because in a 
Markovian bath, the very fast coupling fluctuations suppress the disorder-localization 
effects that hamper the transfer fidelity as we show below for a typical case. 
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Figure 5. (Color online) Transfer infidelity as a function of the transfer time T for 
the optimal control solution in the case of a Markovian bath (Eq. ( [28] ) ) and without 
modulation (p = 0). N + 2 = 31, J\ = Jn- The corresponding boundary-couplings for 
T = 860 are shown in the inset. 


4-3. Non-Markovian noise 


We now consider a non-Markovian noise of the form ./,- + JjAj(f), where A*(£) = 
Aj ([£/t c ]), where the integer part [t/r c ] = n defines a noise A ?; (n) that randomly varies 
between the interval [—£j, £j} at time-intervals of r c during the transfer. We observe 
a convergence of the transfer fidelity to its value without noise as the noise correlation 
time r c decreases (Fig. HR Consequently the fidelity can be substantially improved 
by reducing olm- The effective noise strength scales down as (Fig. ijy inset). 
By contrast to the Markovian limit r c —> 0, dynamical control can strongly reduce the 
infidelity in the non-Markovian regime that lies between the static and Markovian limits 
and whose bath-spectrum is gapped. 


4-4- Other sources of noise 

Timing errors : In addition to resilience to noise affecting the spin-spin couplings, there 
is another important characteristic of the transfer robustness, namely, the length of 
the time window in which high fidelity is obtained. The fidelity F(t) under optimal 
dynamical control (p = 1,2), yields a wider time-window around T where the fidelity 
remains high compared with its static (p = 0) counterpart. This allows more time for 
determining the transferred state or using it for further processing. Consequently, the 
robustness against timing imperfections [551 [53] is increased under optimal dynamical 
control. 

On-site energy noise : This kind of noise, alias diagonal-noise, can be either static or 
fluctuating. The static one can give rise to the emergence of quasi-degenerate central 
states. Then, the dynamical control approach introduced in this work is still capable of 
isolating the “system” defined here (Sec. [TJ from the remaining “bath” levels. It may 
happen that the spin network is not symmetric with respect to the source and target 
spins, and then the effective couplings of the source and target qubits with the central 
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level will not be symmetric. This asymmetry can be effectively eliminated by boundary 
control. On the other hand, a fluctuating diagonal-noise that may produce a fluctuation 
of the central energy level is here fought by optimizing the tradeoff between speed and 
fidelity as detailed above. Additional dynamical control of only the source and target 
spins can be applied to avoid these decoherence effects, by the mapping to an effective 
3-level system, as a variant of dynamical decoupling (56j EH |58) 159| . 

5. Conclusions 

We have proposed a general, optimal dynamical control of the tradeoff between the speed 
and fidelity of qubit-state transfer through the central-energy global mode of a quantum 
channel in the presence of either static or fluctuating noise. Dynamical boundary- 
control has been used to design an optimal spectral filter realizable by universal, 
simple, modulation shapes. The resulting transfer infidelity and/or transfer time can be 
reduced by orders of magnitude, while their robustness against noise on the spin-spin 
couplings is maintained or even improved. Transfer-speed maximization is particularly 
important in our strive to reduce the random phase accumulated during the transfer 
when energy fluctuations (diagonal noise) affect the spins [7BJ. We have shown that, 
counterintuitively, static noise is more detrimental than fluctuating noise on the spin- 
spin couplings. This general approach is applicable to quantum channels that can be 
mapped to Hamiltonians quadratic in bosonic or fermionic operators mmmmm 
We note that our control is complementary to the recently suggested control aimed 
at balancing possible asymmetric detunings of the boundary qubits from the channel 
resonance mm 
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Appendix A. The Hamiltonian in the interaction picture 

The system-bath Hamiltonian (Eq. ([7]) of the main text) splits into a sum of 
symmetric and antisymmetric system operators that are coupled to odd- and even- 

4 

bath modes: H SB (t) = JfSj ® R], where 5/(3) = o(t)(c 0 + (-)c/v+i), *$2(4) = 5{ (3) , 

j = 1 

i?i(3) = Jkbk and £>2(4) = -S|( 3 ). In the interaction picture Hs B {t ) becomes 

^^■^odd(even) 

4 

3 =1 
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where 

Sj{t) = Ul(t)Sj(t)U s (t ), U s (t) = 

Bj(t) = Ug(t)BjU B (t), U B (t) = e~ iHBt - 
and the evolution operators are 

U s (t) = |0) ss (0| + ( c ° s( ^f))+i ) (|0)(0| + \N + 1)(N + 1|) 

+ ( — (|0)(iV + 1| + \N + 1)(0|) 

+ cos( v / 20(t))k)(^| - i sin( ^ (t)) (|0)(z| + | N+ l){z\ + h.c .), 

N 

U B (t)= E e-^|fc)(fc| + |0 )bb<0 |, 

k=l ,k^z 


(A.2) 


(A.3) 


where the states |0)s = |0oCh0jv-i-i)s and |0 )b = | 0 i... 0 jv)b refer to the zero-excitation 
states in the system (S) and bath (B) respectively. Therefore, the bath operators are 

= E iy| 2 e-“*'l*)B{0|, B m (t) = b| w («). 

^^^odd(even) 

We define a basis of operators 0 l to describe the rotating system operators Sj(t) 
via a rotation-matrix 0^(2). They are given by 

— |0)s ((0| + (N + 1|) z> 2 = i>\, 

z>3 = |0)s(^| h = hi, (A.4) 

h = |0) s ((0| - (N + 1|) h 6 = v\, 


such that Sj ( t ) = ' YfOlj ii ( t ) v i . Given that S \( t ) = < fi ( t ) (cos(y / 20(f))hi — i v / 2sin(v / 20(f))z> 3 ), 
2—1 

S 3 ( t ) = 4>( t ) u 5, S 2 ( 4 )( t ) = S \, 3 Jt ) the rotation-matrix vectors are 

= fit ) ( cos ( V 2< f >( t )), 0, —i V 2 sin ( V 2( p ( t )), 0, 0,0) 

= 0(t) (0, cos ( y /2(/>( t )), 0, i \/2sm(v / 20(f)), 0, 0) , A ^ 

Q 3 ,i(£) =0(t)(0,0,0,0,1,0) 1 ' j 

n ^{ t ) = 0(t)(O, 0,0,0,0,1). 


Appendix B. The fidelity in the interaction picture 


Here we derive Eqs. (12 13) from Eq. @ of the main text. Considering \ if ) = 
|100...0)sb = ||0)b with | fj)s = |1o0 z 0jv+i)s as the initial state, the fidelity is 


reduced to 


(B.l) 


fo,N+i(T) = \s (1>\Ps(T ) \1>) s \ = 1 - con , 

where COO = T it > with 

f)g = Si t2 diy + Si )2 8 §y + Sifidiy + (B.2) 


i,i'=l 

r.-.-' = sty I [c,w l , 0)ss 


2,2 


and 


r-T r t 


Ri,i'(T)=- I dt I + ( b . 3 ) 


'0 J 0 








Optimized dynamical control of state transfer through noisy spin chains 


16 


Here by and Oj.j are as defined in Eqs. (A.4 A.5[ ), while the correlation functions are 
*«■(*-*')= E + E (B.4) 

k£k 0 dd k£k even 


This leads to the infidelity £(T) of Eq. (13). 


Appendix C. Considerations for a specific non-Markovian bath: the 
uniform spin-channel 


Consider a uniform (homogeneous) spin-chain channel, i.e. Ji = J in Eq.Q, whose 
energy eigenvalues are 0Uk = 2 Jcos(-M^). In the weak-coupling regime where ajf -C 1, 

the coupling strength in the interaction H & c , J z = and = J 2 sm(^), are 

always much smaller than the nearest eigenvalue gap \u z —oj z ± i| ~ ^ [38, 39, EJ- The 
correlation function of the bath is 


$ ± (r) 


E 

^odd(even) 


7V+1 


Jsin( 


kn 

N+l 


g-*2Jcos(^ T )r 


(C.l) 


and has recurrences and time fluctuations due to mesoscopic revivals, while at short 
times t, it behaves as a Bessel function $(i) = 2 GoP J 1 (2Jt). The latter correlation 
function represents the limiting case of an infinite channel and it gives a continuous 
bath-spectrum that becomes a semicircle. In the case of a finite channel, G(co) will be 
discrete but modulated by the semicircle with a central gap. If disorder is considered, the 
position of the spectrum lines fluctuates from channel to channel but they are essentially 
modulated by the semicircle with a central gap as was considered in the Fig. El 3 of the 
main text, where 

G±{u>) = ^V4J 2 — cu 2 (l - 0(w - (n z )0((n + u> t )), ug = (C.2) 

This is the Wigner-distribution for fully randomized channels [66] with a central gap. 
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